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I.  INTRODUCTION 


One  of  the  fundamental  problems  in  electromagnetic  theory  is  the 
transmission  of  a  spherical  wave  through  a  dielectric  shell.  This  problem 
has  numerous  applications  in  antenna  radomes,  electromagnetic  shielding,  and 
scattering.  It  appears  that  solutions  to  this  problem  are  available  only 
for  the  special  case  where  the  shell  is  an  infinite  dielectric  half  space. 
That  case  was  first  studied  by  Sommerfeld  in  1909,  whereas  later  research 
was  summarized  in  a  book  by  Brekhovskikh  (Chapter  IV  of  [1]).  In  the 
present  paper,  we  consider  a  more  general  case,  namely,  the  shell  has  two 
spherical  boundary  surfaces.  Unlike  the  Sommerfeld 's  problem,  our  case 
does  include  the  effects  of  the  shell's  curvature  and  thickness.  Therefore, 
its  solution  should  be  of  more  practical  interest. 

To  solve  our  problem  rigorously,  the  spherical  wave  expansion  may  be 
used.  However,  due  to  the  fact  that  the  source  location  and  the  two 
dielectric  surface  centers  do  not  coincide,  the  translational  addition 
theorem  for  vector  spherical  wave  function  [2]  must  be  used.  (Our  problem 
is  roughly  comparable  to  scattering  by  three  dielectric  spheres.)  This 
theorem  leads  to  a  complex  series,  which  makes  it  very  difficult  to  generate 
numerical  results.  In  this  paper,  we  use  the  geometrical  optics  theory  (GO) 
[3],  [4]  to  calculate  the  transmitted  field  in  the  problem  sketched  in 
Figure  1.  Such  a  solution,  though  only  approximately  valid  for  high 
frequencies,  is  given  in  a  simple  closed  form.  Thus,  it  allows  us 
to  study  the  "cause  and  effect"  of  the  various  parameters  in  a  convenient 
manner  and  gain  physical  insight. 


II.  SOLUTION  FOR  CONCENTRIC  SPHERICAL  SHELL 


Let  us  consider  the  wave  transmission  problem  in  Figure  1.  The 
boundary  surfaces  of  the  dielectric  shell  are  two  concentric  spheres 
with  radii  R^  and  (both  positive),  and  with  a  common  center  at  Q. 

The  point  source  at  0  emits  a  spherical  wave  described  by  (for  exp  jwt 
time  convention) 

"jkor 

E1  (r)  =  —  [P(e, j)9  +  Q(8,$)»)  .  (2.1) 

Here,  (r,6,<J>)  are  spherical  coordinates  of  r  with  origin  at  0.  Functions 

P  and  Q  describe  the  radiation  pattern  of  the  source.  The  wave  number  . 

1/2 

kg  =  2i r/Ag  =  u)(ygEg)  is  that  of  free  space.  By  using  GO  [3],  [4], 
the  transmitted  field  at  an  observation  point  3  is  to  be  determined. 

Without  loss  of  generality,  we  assume  that  point  3  lies  in  the  (x,z)-plane 
with  rectangular  coordinates  (x^,0,z^). 

A.  Ray  Tracing.  In  accordance  with  Snell's  law,  we  trace  a  ray 
from  source  point  0  to  observation  point  3,  via  refraction  points  1  and 
2  (Figure  1).  Clearly,  all  four  points,  0  to  3,  lie  in  the  same 
(x,z)-plane.  For  a  given  launching  angle  8  of  the  ray  and  the  distance 
c,  the  other  geometrical  parameters  can  be  determined  from  the  following 
relations : 

sin  =  (d/R^)sin  8  ,  a  =  R^[sin(8  -  a*)]/sin  8  (2.2) 

,  t  -1  .  i  .  i  t 

sin  =  n  sm  ,  sm  =  (R^/R2)  sm 

b  ■  R2[sin(a^  -  a^l/sin  ,  sin  =  n  sin 

x^  =*  a  sin  8  +  b  sin(8  -  a^  +  a^)  +  c  sin  (a^  +  +  9) 

3 


'  ■*  _ _ _ _  -  ... 


=  a  cos  0  +  b  cos  (0  -  +  a^)  +  c  cos  (cr  -  a*  -  a*  +  9) 


Thus,  for  a  given  (d,c),  we  can  determine  the  position  of  point  3 
straightforwardly.  On  the  other  hand,  when  point  3  is  given,  explicit 
formulas  do  not  exist  for  determining  0.  One  has  to  find  6  by  trial- 
and-error . 

B.  Field  on  the  Ray.  The  present  vector  field  problem  can  be 
decomposed  into  two  scalar  ones:  one  with  the  electric  field  vector 
perpendicular  the  plane  of  incidence  (E  =  y  E  )  and  the  other  with 
the  electric  field  vector  parallel  to  the  plane  of  incidence  (E  ) . 

-»-t 

The  final  solution  for  the  transmitted  field  E  at  point  3  derived  by 
GO  is  gi^en  by 


EC(3)  = 


E l  (3) 

Tj_  E  l  (1) 

=  (DF) 

=11  (3) 

T„  Ef,  (1) 

L  J 

-jkQ(nb+c) 


(2.3) 


The  various  factors  in  (2.3)  are  explained  below:  Tj_  and  are  the 
products  of  the  transmission  coefficients  at  points  1  and  2,  given  by 
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(2.4b) 


The  two  components  of  the  incident  field  are  calculated  from  ( 2 . 1) : 

-jkoa 

EiL(l)  =  - -  Q(e,$  =  0)  (2.5a) 


-jk0a 


E,rd)  = 


P(e,4>  =  o) 


(2.5b) 


-V  t 


The  derivation  of  the  divergence  factor  DF  in  (2.3)  requires  some  effort. 
With  the  details  given  in  Appendix  A,  the  final  result  for  DF  reads 


DF  =  (1  +  <11b)"1/2  (1  +  K12b)"1/2  +  ,c21c)~1/2  (1  +  k22c)"1/2  (2*6) 


At  point  1  or  2,  it  can  be  shown  that  the  principal  directions  of  the 
transmitted  wavefront  are  precisely  the  two  directions  parallel  and 
perpendicular  to  the  plane  of  incidence.  Here  (K21,IC12^  are  tbe  tw0 
principal  curvatures  of  the  transmitted  wavefront  at  point  1,  and 
(<2i>K22^  are  tbose  at  Point  2.  They  are  calculated  from  the  relations 


=  (n  cos  a 


t.-lfl  2  i  1  ,  t 

)  —  cos  a  +  — (n  cos  ctn 

13.  1  1 


1,1,  t  1  iN 

c10  =  —  +  -r—  (cos  a, - cos  a,) 

12  na  R^  In  1 


(2.7a) 

(2.7b) 


,  2  t,-l  ..  -l.-l  2  i  .  1  ,  t  i., 

«2i  “  (cos  a2)  I  (b  +  7  n  cos  +  ^(cos  a2  -  n  cos  a2)J 


(2.7 c) 


k22  =  n(b  +  <12)  1  +  l-(cos  ab  -  n  cos  a2)  .  (2.7d) 

The  sign  convention  of  k  is  as  follows.  If  ic  is  positive  (negative), 
the  normal  section  of  the  wavefront  is  divergent  (convergent) .  For 


example,  if  the  transmitted  wavefront  at  point  1  is  the  same  as  the 

incident  spherical  wavefront,  we  have  k  ^  =  k^2  =  +R^.  For  a  typical 

factor  in  (2.6),  the  square  root  convention  is 

+| f { ,  if  f  is  real 

f  =  (1  +  Kb)'1/2  =  ■  •  (2.8) 

+ j  | f |  ,  if  f  is  imaginary 


When  f  is  imaginary,  it  means  that  the  ray  has  crossed  a  focus  of  the 
ray  pencil.  The  (+j)  accounts  for  the  well-known  (n/2)  phase  retardation. 


Several  general  comments  about  the  solution  in  (2.3)  are  in  order: 

(i)  For  the  case  where  a  total  reflection  occurs  (a^  or  becomes  complex), 
the  field  in  the  transmitted  region  is  not  an  optical  field,  and  the 
present  ray  solution  (2.3)  is  no  longer  valid,  (ii)  Except  for  special 
cases,  e.g.,  normal  incidence  =  0,  the  two  curvatures  (,c2i,,e22^ 
the  transmitted  wavefront  emerging  from  the  dielectric  shell  are  not 
equal.  Thus,  the  transmitted  pencil  is  generally  astigmatic.  (iii)  It 
is  possible  that  and/or  are  negative.  Then  the  divergent  pencil 
from  the  source  is  transformed  into  a  convergent  (focusing)  pencil  after 
propagating  through  the  dielectric  shell.  (iv)  The  solution  in  (2.3) 
remains  valid  for  more  geometries  than  the  one  shown  in  Figure  1.  This 
is  discussed  further  in  the  next  section. 

III.  GENERALIZATION  OF  FINAL  SOLUTION 

For  transmission  through  the  concentric  spherical  shell  in  Figure  1, 
the  final  solution  consists  of  two  parts: 

Part  A:  Ray  tracing  formulas  in  (2.2) 

Part  B:  Field  solution  in  (2.3)  through  (2.5) 

It  can  be  shown  that  Part  B  is  valid  under  a  more  general  condition 
(Figure  2),  namely, 

(i)  The  centers  of  the  spherical  shells  Q1  and  Q?  need  not 

coincide,  as  long  as  the  four  points  (0,1, 2, 3)  are  coplanar. 

(ii)  The  surfaces  of  the  shell  can  be  either  concave  or  convex. 

Looking  from  the  source  side,  R^  (or  R2)  is  positive  if  the 
surface  is  concave,  and  is  negative  if  the  surface  is 
convex. 


.  v 
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Condition  (i)  is  necessary  in  order  for  the  scalar  transmission  coefficients 
and  Tj (  in  (2.4)  to  be  valid.  When  the  four  points  (0, 1,2,3)  are  not 
coplanar,  the  perpendicular  and  parallel  polarizations  are  no  longer 
uncoupled.  Then  the  scalar  Tj_  or  T 1 1  in  (2.3)  must  be  replaced  by  a 
(2  x  2)  matrix.  Also,  (2.7)  becomes  more  complicated. 

For  the  general  configurations  in  Figure  2,  the  ray  tracing 
formulas  in  (2.2)  are  not  valid.  However,  by  following  Snell's  law, 
the  ray  tracing  even  in  the  most  general  situation  is  conceptually  simple. 
Thus,  instead  of  working  out  a  set  of  general  formulas,  we  leave  it  to  the 
individual  problems. 

IV.  AXIAL  INCIDENCE  ON  SYMMETRICAL  SHELLS 


To  study  the  features  of  the  present  ray  solution,  let  us  concentrate 
on  a  special  case,  where  the  four  points  (Q^.C^.O^)  are  on  a  straight 
line  (Figure  3).  Then  the  four  curvatures  in  (2.7)  reduce  to 


1  ,  n  -  1 

<,  ,  =  <,  »  =  -  +  - — 

11  12  na  n  R, 


<22  *  <21 


nR^  +  an(n  -  1) 
bR^  +  ab(n  -  1)  +  na  R^ 


1  -  n 


n 


R~  .  ,  ~1  R~ 

2  b  +  <u  2 


(4.1a) 


+  1  I  n  >  (4.1b) 


An  interesting  question  is  "When  does  <^  or  <22  become  negative  (meaning 
a  convergent  circle)?"  This  is  answered  below: 

(i)  Negative  <^.  The  Cransmitted  pencil  inside  the  dielectric 


shell  is  a  convergent  one  when  <^  <  0  or 
1 


a  > 


1  -  n  “1 


>  0 


(4.2) 


If  R^  >  0  (concave  dielectric  interface  shown  in  Figure  4a),  this  is  possible 
if  n  <  1.  If  <  0  (convex  dielectric  interface  shown  in  Figure  4b), 


earn  inside  the 


this  is  possible  only  if  n  >  1.  The  distance  between  point  1  and  focal 
point  is  (k^) 

(ii)  Negative  <22*  The  transmitted  field  in  the  free-space  region 
outside  the  dielectric  shell  is  convergent  when  <22  <  Let  us  concentrate 
on  a  special  case  in  which  the  thickness  of  the  dielectric  shell  is  small 
so  that 


b  <<  k 


-1 

11 


(4.3) 


Then  < ^  in  (4.1b)  becomes  approximately 


<22  :  7  +  (n  -  1) 


tRl 


1_ 

R„ 


(4.4) 


which  is  the  well-known  equation  for  the  thin  lens.  (See  for  example 
Eq.  (41-1),  p.  685  of  [5].  Note  the  corresponding  notations  used  in  [5] 
and  here:  s  *  -a,  s'  -  *  ~Ri>  an<*  Ro  “  ~R->-)  The  condition 

for  a  negative  <22  under  the  approximation  in  (4.3)  is 


a  > 


R1R2 


(n  -  1)0^  -  R2) 


>  0 


(4.5) 


(iii)  Far  Field:  If  the  observation  point  3  is  in  the  far  zone 


(<29c  >>  1),  then  DF  in  (2.6)  becomes 


DF 


a  +  —  b  +  c 
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1  +  (n  -  1) 


ab 

R1R2 


R~  -  R, 


1  +  - 
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i-i 
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h(4.6) 


The  first  factor  in  {  }  in  (4.6)  is  the  divergence  factor  of  a  planar 
dielectric  slab  (R^,R2  —  ®) .  Thus,  the  ratio  of  the  electric  field  at 
a  far-fleld  point  3  for  a  spherical  shell  and  that  for  a  dielectric 
slab  is 


11 
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|1EC(3)  for  spherical  shell] _ 

|eC(3)  for  a  slab  of  same  thickness) 


1  +  (n  -  1) 


ab 


rir2L 


R2  -  ri  i 

— =—r — -  -  1  +  - 
b  n 


i-i 
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(4.7) 


As  a  numerical  example,  consider  the  case  in  which  the  inner  and  outer 
dielectric  surfaces  are  concave  and  concentric  (Q^  =  Q^)  with  (R^/b)  «  2. 
We  plot  n  as  a  function  of  (a/b)  for  n  *  0.5  and  n  *  3  in  Figure  5.  We 
note  that  n  can  be  substantially  different  from  unity.  When  a  =  R^,  we 
have  r\  *  1.  Thus,  in  this  interesting  special  case,  the  axial  far  field 
through  a  concentric  spherical  shell  and  that  through  a  planar  slab 
become  the  same.  Another  interesting  special  case  occurs  when  DF  -*•  °°. 

It  means  that  the  paraxial  rays  emerging  from  the  dielectric  shell 
(Figure  3)  are  parallel  to  the  axis  so  that  they  focus  at  the  far-field 


point  at  infinity.  From  (4.6),  it  is  shown  that  DF  -  °°  if  < 

^  R^  -  R2  +  b(n  -  l)/n 

a  *  (n  “  X)  R1[R2  -  b (n  -  l)/n] 


22 


0  or 


(4.8) 


Under  the  thin-lens  approximation  b  ;  0,  (4.8)  is  reduced  to  the  well- 
known  lensmaker's  equation  (see  Eq.  (41-2),  p.  685  of  [5]).  In  the 
antenna  radome  application,  (4.8)  is  useful  in  the  determination  of  the 
enhancement  of  the  antenna  main  beam. 

(iv)  Multiple  Refraction:  For  a  given  source  point  0  and  observation 
point  3  in  Figure  1  or  2,  we  can  trace  two  types  of  geometrical  optics  rays 
The  first  type  is  the  direct  ray  from  0  to  3  without  going  through  internal 
reflections  in  the  dielectric  shell.  Its  field  solution  is  given  in  (2.3) 
which,  of  course,  is  the  main  contribution.  The  second  type  contains 
rays  which  bounce  one  or  more  times  inside  the  shell  before  reaching  point 
We  now  consider  the  contribution  of  such  multiply  refracted  ravs.  For  the 


axial  incidence  case  (Figure  3),  let  us  further  specialize  the  geometry 


slab 


1  - 


n  -  1 


exp(-j2kQnb) 


(4.12c) 


ln  +  1J 

Comparing  (4.12b)  with  (4.9),  we  note  the  effect  of  the  multiply  refracted 
rays  is  accounted  for  by  replacing  T  by  T  ^en  the  condition  in 

(4.12a)  is  not  met,  we  must  evaluate  (4.11)  numerically.  Let  us  define 


an  error  term 


6  = 
P 


|E  (3) |  including  p  multiply  refracted  rays 


-  1 


idrt 


[|E  (3) |  including  all  multiply  refracted  rays 


x  100%  .  (4.13) 


In  Figure  7,  we  plot  6^  (including  no  multiply  refracted  rays)  and  6^ 
(including  one  multiply  refracted  ray)  vs.  n  for  b  =  0.75  X^  and 
(a  +  c)/b  =  4.  Several  observations  are  made.  (a)  For  commonly  used 
values  of  n  (between  1  and  3) ,  the  error  <$q  is  13%  or  less  except  at 
resonances.  A  resonance  occurs  when  all  multiply  refracted  rays  emerging 
from  the  slab  are  in  phase  with  the  primary  ray.  For  the  configuration 
in  Figure  6,  the  resonance  condition  is 


(nk^b/ir)  =  (2nb/XQ)  =  a  positive  integer 


(4.14) 


For  an  obliquely  incident  ray  and/or  a  curved  slab,  the  condition  for 
resonance  is  rarely  satisfied.  Thus,  generally  speaking,  the  error  for 
neglecting  the  multiple  refraction  is  roughly  10%.  (b)  Errors  6^  and  6^ 

have  about  the  same  order  of  magnitude.  Thus,  the  inclusion  of  the  first- 
order  multiply  refracted  ray  does  not  in  general  improve  the  accuracv  of 
the  solution. 
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V.  NUMERICAL  RESULTS 


A.  Radome .  Consider  a  dielectric  spherical  radome  with  an  inner 
radius  20  Aq,  a  uniform  thickness  0.5  A^  and  n  =  /5",  as  shown  in  Figure  8. 
The  y-polarized  source  is  located  at  points  1,  3,  or  4,  and  its  radiation 
field  in  the  E-plane  is  given  by 

-iR0r 


EX  (r ,  0  ,  <J)  =  tt/2)  =  - — - —  [y  cos  (1.5  6)] 


(5.1) 


The  transmitted  field  is  calculated  from  (2.3).  For  an  observation  point 
in  the  E-plane  and  in  the  far-field  2one  (x  =  0  plane  and  c  ■»  in  Figure  1) , 
we  may  rewrite  (2.3)  as 

"jV  „ 

|c(r)  '  - — - —  [y  pt(6)]  ,  r  -*■  “  (5.2) 

where  r  is  the  distance  from  point  0  to  point  3  (Figure  1).  We  plot  pC(6) 
as  a  function  of  8  in  Figure  8.  Generally  speaking,  the  radome  modifies 
the  radiation  field  gently,  as  expected. 

B.  Lens  I  (Double  Concave).  Unlike  the  above  radome,  a  dielectric 
lens  may  modify  the  incident  field  drastically.  Let  us  consider  Lens  I, 
drawn  approximately  to  scale  in  Figure  9.  The  source  is  2  \  away  from 
the  lens,  and  is  y-polarized.  Let  us  concentrate  on  the  field  in  the 
H-plane  (x-z  plane).  In  Figure  9,  we  launch  4  rays  2°  apart.  The  outside 
ray  (at  9  =  6°)  suffers  total  reflection  at  the  second  face  of  the  lens, 
and  is  not  transmitted  into  the  free  space  region  (we  ignore  multiple 
refractions).  We  assume  the  incident  field  from  the  source  is  confined 
to  a  cone  (a  beam).  In  the  x-z  plane,  it  is  given  by 


Ei(x,0,z) 


if  9  <  e1 
if  e  >  e1 


(5. 3") 
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In  the  present  case  (Lens  I),  we  choose  0*  =  5.5°.  After  transmission 
through  the  lens,  the  rays  become  more  divergent,  and  they  are  spread 
over  an  angular  region  of  about  0  <  53.5°,  which  is  considerably  wider 
than  the  incident  angular  spread  (0  <  5.5°).  At  an  observation  point 
(r,6,<)>  *  0),  we  express  the  transmitted  field  as 

_jkor  . 

EC(r,0,O)  =  -e- . —  [y  QC (r , 0) ]  .  (5.4) 

In  the  absence  of  the  lens,  QC  =  1  for  0  <  S*  and  QC  =  0  for  0  >  0*. 

t  3 

With  the  lens  present,  we  plot  Q  (r,0)  as  a  function  of  0  for  r  =  2  x  10  X^ 

in  Figure  10.  Note  that  the  transmitted  field  is  much  weaker  (12%  or  less) 

than  the  incident  field  because  of  the  wider  spread  of  the  transmitted  rays. 

For  the  present  case,  QC(r,0)  is  only  very  weakly  dependent  on  r,  as 

long  as  r  >  100  X^.  Thus,  the  transmitted  field  E*"  in  (5.4)  in  the  far 

zone  is  approximately  a  spherical  wave  with  an  angular  pattern  QC.  Every 

transmitted  ray  has  two  foci.  Their  distance  behind  the  second  face 

of  the  lens  is  (k^)  an(*  (*22^  ’  may  be  calculated  from  (2.7). 

In  particular,  (k^)  ^  is  f°r  t'ie  n°rmal  section  of  the  wavefront  in  the 

plane  of  incidence  (x-z  plane),  whereas  (k^)  ^  is  ^or  c^at  in  the 

perpendicular  plane  (defined  by  the  y-axis  and  the  ray  direction) .  For 

the  third  ray  (incident  0=4°)  in  Figure  9,  we  calculate  from  (2.7)  that 

(k21)"L  "  +0-21  *0  »  (|C22)_1  =  +0'375  V  (5'5) 

We  mark  the  position  of  the  focus  corresponding  to  by  a  cross  in 
Figure  9.  The  distance  between  A  and  the  cross  is  (<2^  The  trace 
of  the  two  sets  of  foci  is  shown  in  Figure  11.  They  are  curves  on  the 
two  caustic  surfaces  of  the  transmitted  wavefront  (intersection  of  caustic 
surfaces  and  the  x-z  plane). 
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C.  Lens  II  (Double  Convex).  The  geometry  of  Lens  II  and  its  ray 
picture  are  shown  in  Figure  12.  The  incident  field  is  given  in  (5.3)  with 
9*  *  10°.  The  transmitted  field  in  the  H-plane  is  expressed  in  (5.4),  where 
Qt(r,0)  is  again  very  weakly  dependent  on  r  in  the  far  zone  and  is  plotted 
in  Figure  10.  We  note  that  the  incident  10°-beam  is  now  spread  into  a 
67°-beam  after  transmission  through  Lens  II.  The  peak  value  of  the 
transmitted  field  is  about  242  of  the  incident  field  at  the  same  far-field 
location.  The  caustic  curves  are  shown  in  Figure  11. 

D.  Lens  III  (Convexo-Concave) .  The  geometry  of  Lens  III  is  shown 
in  Figure  13.  The  radii  (R^.R^)  of  the  lens  and  the  source  distance 
satisfy  the  lensmaker's  equation  in  (4.8),  so  that  the  transmitted  rays 
near  the  axial  direction  are  almost  parallel  and  focus  at  a  point  at 
infinity  in  the  axial  direction.  The  incident  field  is  given  in  (5.3) 
with  01  *  12°.  The  H-plane  transmitted  field  E1"  is  expressed  in  (5.4), 
where  Qt(r,0)  varies  drastically  from  the  near  field  zone  up  to  r  '  100  \ 
as  seen  from  Figure  14.  Beyond  r  '  100  the  beam  becomes  narrower,  and 
the  peak  becomes  higher  with  the  increase  of  r.  It  is  well-known  that 
the  exact  value  of  the  beam's  peak  (on  a  caustic  surface)  cannot  be 
predicted  by  the  present  geometrical  optics  theory.  It  can  be  calculated 
from,  for  e.g.,  the  Huygens-Green  formula  described  in  p.  107  of  [7]. 

The  caustic  curves  are  shown  in  Figure  15. 
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Crosses  indicate  foci. 


which  ] ie  in 


VI.  CONCLUSION 


For  a  given  incident  field  in  (2.1),  a  geometrical  optics  solution  of 

the  transmitted  field  through  a  spherical  dielectric  shell  (Figure  2)  is 

given  in  (2.3)  through  (2.7).  The  present  solution  is  approximately  valid 

when  (i)  the  source  frequency  is  high  (the  radii  of  curvature  of  the  curved 

* 

shell  are  large  in  terms  of  wavelength),  (ii)  the  four  points  (0,1, 2, 3)  in 
Figure  2  are  coplanar,  and  (iii)  the  observation  point  is  away  from  the 
caustic  surfaces  of  the  transmitted  wavefront.  Furthermore,  our  solution 
includes  only  the  contribution  from  the  first-order  refracted  rav  (a  direct 
ray  from  '‘he  source  to  the  observation  point  without  internal  refractions 
inside  the  dielectric  shell).  It  has  been  estimated  that,  except  at 
"resonances"  which  rarely  occur  for  curved  shells,  the  error  of  neglecting 
higher-order  refracted  rays  is  roughly  10T.  for  n  =  in  the  range  1  to  3. 
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APPENDIX  A 


DERIVATION  OF  EQUATION  (2.6) 


Referring  to  Figure  1,  the  curvature  matrix  [3],  [6]  of  the  incident 
spherical  wave  may  be  expressed  as 
-1 

U  I 

(A.  1) 


Qi  = 


-l 


The  curvature  matrix  of  the  inner  surface  of  the  radome  (£^)  is 


0 


(A. 2) 


From  the  knowledge  of  these  two  curvature  matrices,  the  curvature  matrix  of 
the  refracted  ray  1-2  may  be  expressed  as  [3] 

I, 


^i  =  ^ei^  +  (cos  ai  ~  cos  ai/n)Q  1Hp^)  1 

where  8^  and  9^  are  coordinate  transformation  matrices  given  by 


cos  a. 


(A.  3) 


(A. 4) 


*1  " 


cos  a. 


(A. 5) 


Simplifying  (A. 3),  we  obtain 
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(A. 6) 


where  and  are  defined  in  (2.7). 

The  curvature  matrix  of  the  wavefront  incident  at  point  2  is  given  by 


Q2  =  kqJ)  1  +  b  n  1 


(A. 7) 


where  I  is  the  identity  matrix.  The  curvature  matrix  of  the  transmitted 
wavefront  at  2  can  be  calculated  similarly  to  that  at  1.  The  final  result 


(A. 3) 


where  k ^  and  k99  are  defined  in  (2.7). 
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